function g=find_g_pun(x)
c=[0.201232428 614230.043 996439.097 3.33560562 25.2096377 -0.381360867];
c1=c(1);
c2=c(2);
c3=c(3);
c4=c(4);
c5=c(5);
h=c(6);
x2=(x-h).^2;
g=c1+c2*x2./(c3+x2).*(1+c4*exp(-c5*x2));